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ABSTRACT 

The evolution of the magnetic field in an accreting neutron star is investigated using a 
fully general relativistic treatment and assuming that initially the currents supporting 
the field are completely confined to the crust. We find that the field decay slows down 
due to the inclusion of the curvature of space-time but the final results do not differ 
significantly from those obtained assuming a flat space-time. We also find that such 
modifications are small compared to the uncertainties introduced by a lack of precise 
knowledge of the neutron star micro-physics. 
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1 INTRODUCTION 



The generation and evolution of the magnetic field in neu- 
tron stars continue to evoke much interest, particularly be- 
cause there is as yet no consensus on the question of field 
generation. The field could be a fossil remnant from the pro- 
genitor star in the form of Abrikosov fluxoids embedded in 
the proton superconductor of the core (Baym, Pethick & 
Pines 1969, Ruderman 1972). Or it could be generated after 
the formation of a neutron star in which case the currents 
would be entirely confined to the solid crust (Blandford, 
Applegate & Hernquist 1983, Urpin, Levshakov & Yakovlev 
1986). Alternatively, a seed magnetic field could be ampli- 
fied due to the entropy-driven convection in a young neutron 
star (Thompson & Duncan 1993). Observations and statis- 
tical analyses of existing pulsar data, however, indicate that 
significant decay of magnetic field is achieved only if the 
neutron star is a member of an interacting binary (see Bhat- 
tacharya 1995, 1996, van den Heuvel 1995, Verbunt & van 
den Heuvel 1995 and references therein) irrespective of the 
generation mechanism. 

The processes responsible for the field evolution, in neutron 
stars which are members of binaries, are - a) expulsion of 
the magnetic flux from the superconducting core during the 
phase of propeller spin-down, b) rapid ohmic decay of the 
crustal field in an accretion-heated crust and c) screening 
of the field by accreted matter. Amongst these, the diamag- 
netic screening of the field by accreted matter depends on 
the detailed geometry of the material flow and whether it 
can have a long-term effect is still not very clear (Konar 
1997, Choudhuri & Konar 2002). The rest of the models in- 



voke ohmic decay of the current loops, residing in the crust, 
for a permanent decrease in the field strength (see Konar 
& Bhattacharya 2001 for a brief review). In either case, the 
effect of accretion is two-fold. The heating reduces the elec- 
trical conductivity, and consequently the ohmic decay time- 
scale, inducing a faster decay. At the same time the material 
movement, caused by the deposition of matter on top of the 
crust, pushes the original current carrying layers into deeper 
and denser regions where the higher conductivity slows the 
decay down. 

Though the problem of the evolution of the magnetic field 
in neutron stars have been investigated by many authors a 
fully general relativistic treatment has not been attempted 
by many. Recently, Sengupta (1997, 1998) addressed this is- 
sue for isolated pulsars and found that the modification in- 
troduced by a general relativistic treatment is significantly 
large. However, our results are at variance with this claim 
(Ramachandra & Vishweshwara 2000). Indeed this is sup- 
ported by the recent work of Page, Geppert & Zannias 
(2000) and Geppert, Page & Zannias (2000). These authors 
have investigated coupled spin, magnetic field and thermal 
evolution of isolated neutron stars in a general relativistic 
framework. And they find that the change in the final field 
strength is quite small. 

In a series of papers (Konar & Bhattacharya 1997, 1999a, 
1999b) we have discussed the accretion induced field evolu- 
tion in neutron stars. The formalism adopted in these pa- 
pers implicitly assumes a flat space-time. In this work we 
report the results of our investigation of the field evolution 
in an accreting neutron star adopting a general relativistic 
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framework. Our aim is to calculate the amount of modifica- 
tion introduced by assuming a curved space-time vis-a-vis 
the uncertainties already existing in the field evolution cal- 
culation assuming a flat space-time. It is well known that 
such uncertainties arise out of a lack of knowledge regarding 
the values of various crustal parameters (like the density at 
which the currents are concentrated or the impurity content 
of the crust) and the observations are still not adequate to 
define those with better precision. 

In this work we assume that the currents supporting the field 
are initially completely confined to the crust. This assump- 
tion does not affect the general qualitative conclusions ob- 
tained from our present calculation. Because,it has already 
been shown (Konar & Bhattacharya 1999b) that an initial 
configuration supported by the magnetic fluxoids in the pro- 
ton superconductor behaves in a similar fashion once it is 
converted to crustal currents by spin-down driven flux ex- 
pulsion. 

We also assume the mass accretion to be spherical. Of 
course, except for very large values of mass transfer the 
accretion would be channeled through the polar cap. The 
detailed nature of this material flow and its effect on the 
magnetic field has recently been investigated (Choudhuri & 
Konar 2002, Konar & Choudhuri 2002). But such flows are 
confined only to the layers very close to the surface. But 
the currents would normally be concentrated in the deeper 
curtal layers for the magnetic field to be long lived. And in 
these layers the material movement would be purely radial 
due to the overall compression of the neutron star, as a re- 
sult of an increase in the mass. Therefore, an assumption 
of a spherical accretion would not introduce any significant 
error in our present calculations. 

The organisation of the paper is as follows. In section ^ wc 
derive the covariant form of the induction equation that gov- 
erns the temporal evolution of the magnetic field. In section 
|§] we discuss the detailed form of the induction equation, 
both for a flat and a curved space-time, in the crust of a 
neutron star. Section ^ provides some details of the micro- 
physics of the neutron star required for our investigation. In 
section ^| we present the results of our calculation and finally 
conclude in section []. 



2 THE INDUCTION EQUATION 

The induction equation in a flat space-time is given by (Jack- 
son 1975): 



dB ~ c 2 - 1 

= V x (V x B) - 4-V x(-xVxB). 
dt y ' 4tt V ; 



(1) 



where V is the velocity of material movement and a is the 
electrical conductivity of the medium. To obtain the covari- 
ant form of the induction equation we make use of the co- 
variant form of the Maxwell's equations (Weinberg 1972) : 




4tt 



J", 



dpF"* + duF + d x F 



(2) 
(3) 



with 

TP 

* [If 



(cp, J) 



-d v A, 



where A M is the vector potential, p is the charge density 
and J is the spatial part of the current density. Using the 
generalised Ohm's law given by (Weinberg 1972): 

r = ag^F„ x u X , (4) 

where u A is the covariant velocity, from the first Maxwell's 
equation we obtain: 

c 1 u k 
Fao = ■^gavdviV^F'"') - F ak — . 



(5) 

With the covariant velocity given by u M = |«°V), 
where V is the velocity in the locally inertial frame, we ob- 
tain from the second Maxwell's equation: 

d t Fij = di(F ]k v k ) - dj(F ik v ) 



(6) 




where = (ct,r,8,(f>). It should be noted that we have 
neglected the displacement current here and i, j, k, I = 1, 2, 3 
since we are only interested in the time evolution of the 
magnetic field. For a given F^v an observer with four velocity 
u M measures the electric and the magnetic fields (E,B) as 
given by: 



Efj, = F^ v u and £y 



(7) 



where e M „ 7< 5 is the four-dimensional Levi-Civita tensor ( Wald 
1984). Therefore eq.(^|) gives the covariant form of the in- 
duction equation which reduces to eq.(Q) in the limit of a 
flat metric. We must also mention here that in this work 
we assume that the magnetic field is weak enough for the 
electrical conductivity to be isotropic. 



3 EVOLUTION IN THE CRUST OF A 
NEUTRON STAR 

In this work we shall follow the formalism developed by 
Konar & Bhattacharya (1997 - paper I hereafter) for the 
evolution of a crustal field assuming a flat space-time. And 
shall compare the results with those obtained in paper I to 
obtain the modifications introduced by assuming a curved 
space-time. 

We assume the star to be accreting at an uniform rate, with- 
out any loss of generality, since the purpose of this work is to 
find the modification introduced by the curvature of space- 
time. Therefore we shall, for the moment, ignore the realistic 
situation of a neutron star accreting at different rates as its 
binary companion goes through different stages of evolution. 
Assuming the mass flow to be spherically symmetric in the 
crustal layers of interest, the velocity of the material move- 
ment is given by (see paper I for details): 



v(r) 



M 



4:TY 2 p(r) 



(8) 
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where M is the rate of mass accretion and p(r) is the density 
as a function of radius r. 

The mass in the crust of the neutron star is quite small com- 
pared to the total mass of the star. For example for a inter- 
mediate equation of state (Wiringa, Fiks & Fabrocini 1988 
matched to Negele & Vautherin 1973 and Baym, Pethick & 
Sutherland 1971) the mass of the crust of a 1.4 M^neutron 
star would be 0.044 Mq, which is ~ 3% of the total mass. 
Therefore, in the crust of a non-rotating neutron star we 
can use the exterior Schwarzschild metric given by (Wein- 
berg 1972): 



= l-(i ),(i ) , 

V r r 



without making any significant error. Here, m 
being the mass of the star. 



(9) 



M 



It should be noted here that the crustal mass of a neutron 
star remains effectively constant for the accreted masses 
of the order of 0.1 MQ(in fact, with a slight decrease) as 
the total mass increases. Therefore accretion implies as- 
similation of the original crust into the superconducting 
core. We assume that when the original current carrying 
regions undergo such assimilation the superconducting tran- 
sition happnes over a time-scale much smaller than the flux- 
expulsion time-scale and hence the flux remains frozen there- 
after stopping any further decay. Since we terminate our cal- 
culation long before the crustal mass changes significantly 
(as the field 'levels off' due to flux freezing much before that) 
the use of exterior Schwarzschild metric remains valid even 
in the crust of an accreting neutron star. 



3.1 Flat Space-Time 

In a flat space-time : 
B = V x A, 



(10) 



where we assume A = (0,0, A/,). This choice of the vector 
potential ensures a poloidal geometry for B. In particular, 
we take, 



Aa 



g(r, t) sin 9 



(11) 



where g(r,t) is the Stokes' function. In the lowest order of 
multipole, the dipolar form of the magnetic field is given by: 



B r = 



Bo 



2 cos 9g(r,t) 

sin 9 dg(r, t) 
r dr 



(12) 
(13) 



where, the quantities A, B correspond to those in a flat 
space-time. The induction equation then takes the form: 



dt9(r,t) = — <d r g(r,t) — 



+ d r g(r,t)v(r). (14) 



At the pole the magnitude of the magnetic field is propor- 
tional to 2g ^'^ . Therefore, to find the time-evolution of the 
strength of the surface magnetic field at the pole we solve 
eq.(|14) subject to the following boundary conditions (see 



e.g., Geppert & Urpin 1994): 
dg(r,t) l ,g(R,t) 

=R + 



dr 



R 

g(rco,i) 



(15) 
(16) 



where R is radius of the star and r co is that radius to which 
the original boundary between the core and the crust is 
pushed to, due to accretion, at any point of time. The first 
condition matches the interior field to an external dipole 
configuration. The second condition indicates that as accre- 
tion proceeds along with the crustal material the frozen-in 
flux moves inside the core. For the details regarding the par- 
ticular form of g(r, 0) and the current density corresponding 
to it, please see paper I. 



3.2 Curved Space-Time 

In a Schwarzschild metric F^v is related to the locally mea- 
sured magnetic magnetic field B by the following relations 
(Wasserman & Shapiro 1983): 



(l-2m/r) 1/2 



F, 



B r 



r sm ( 
1 



4>ri 



(17) 
(18) 



Then the choice of, A$ — g(r, t) sin 2 9, ensures that B and B 
match in a locally flat space-time where B is obtained from 
eq.(0). With this choice, we have (Wasserman & Shapiro 
1983): 



Bo 



2 cos 6» 



g(r,t) , 



, 2m \ 1/2 sin(9 n . . 

1 d r g(r,t). 

r ) r 



(19) 
(20) 



Therefore, in a curved space-time the induction equation 
reduces to: 



c t s c 2 / 2m\!/ 2 



1 - — ) &rg(r,t) 



+ 



c 2 / 2m\ 1 / 2 2m 

i — i + 

47TO- V r ) r 2 



d r g(r,t). 

(21) 



Once again the strength of the field at the pole is simply 
proportional to g(R,t). But now the boundary conditions 
get modified to the following form (see Zannias, Geppert & 
Page 2000): 

d r g(r,t)\ r =H + g i ^- = 0, (22) 
ff(roo,t) = 0; (23) 

with, 

2yMl -y- 1 ) + ^ 

Q = y^Tir — ~itt r 1/9 > y = R l 2m - ( 24 ) 

y 2 m(l-y 1 )+y + l/2 



4 NEUTRON STAR MICRO-PHYSICS 
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M/M 


P71 


WFF 2 


W74 


1.4 


7.45 


11.08 


12.21 


1.6 




11.03 


12.18 


1.8 




10.91 


12.10 


2.0 






11.91 



Table 1. Radii of neutron stars of different masses for different 
equations of state, in Km. The equations of state P71, WFF and 
W74 have been discussed in the text. 

4.1 Equation of State 

In order to investigate the effect of the curvature of space- 
time on the evolution a purely crustal field, we consider three 
equations of state - soft, intermediate and stiff. The soft eos 
is taken from Pandharipande (1971) (P71 hereafter), the in- 
termediate from Wiringa, Fiks & Fabrocini (1988) (WFF 
hereafter) and the stiff eos from the work by Walecka (1974) 
(W74 hereafter) respectively. All these eos are matched to 
Negele & Vautherin (1973) and Baym, Pethick & Suther- 
land for the low-density crustal region. In table (Q) we have 
listed the radii of stars of different masses obtained by using 
different equations of state. 
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Figure 1. The evolution of the surface magnetic field in an 
accreting neutron star of mass 1.4 MQCorresponding to an in- 
termediate equation of state given by WFF, assuming a flat 
space-time. Legends next to the curves - 1 and 2 correspond to 
M = 10 — 12 ,10 — 10 Mo/yr respectively, A and B correspond to 
densities of 10 14 ,10 13 gem -3 respectively (at which the initial 
current configurations are centred) and a and b correspond to 
impurity contents of and 0.01 respectively. 



4.2 Crustal Physics 



with, 



As has been discussed in paper I for the uniform crustal 
temperature of a neutron star, accreting at a rate of M, we 
use a fitting formula to the results obtained by Zdunik et al. 
(1992) given by: 



log T = 0.397 log M + 12.35 . 



(25) 



Extrapolation of this fit to higher rates of accretion gives un- 
physically high temperatures. It has been shown by Brown 
(2000) that the maximum temperature obtained in an accre- 
tion heated crust is ~ 10 8 ' 5 K. Therefore we use this value 
if the temperature obtained through eq.(p5|) exceeds it. It 
should be noted here that the above results were obtained 
for a neutron of mass 1.4 MQwith an intermediate equation 
of state. Since, such calculations do not exist for different 
masses or different equations of state, use of eq.(p5|) for all 
situations introduce some error in our calculation. But since 
the purpose of this work is to obtain a measure of the mod- 
ification in the final field strength by using a general rela- 
tivistic framework we shall, for the moment, use eq.(^) as 
an approximate indicator for the crustal temperature. 

In the solid crust, the electrical conductivity has contribu- 
tions from both the phonon and the impurity processes and 
is obtained as, 

i = -U-^ (26) 

O <7 p h CTi mp 

where <r p h is the phonon scattering conductivity dependent 
upon the density and the temperature (Itoh et al. 1984) 
and (Jimp is the impurity conductivity dependent upon the 
density and the impurity parameter Q (Yakovlev & Urpin 
1980). Whereas, a p ^ and o"i mp are given by, 



c p h 



1.24 x irr 



) x 4 (u 2 + Q.Q174) 1 / 2 
uT s (1 + 1.0183; 2 )/,, 
8.53 x W^xZ/Qs' 1 



(27) 
(28) 



Ts = temperature in units of 10 8 K 
pe = density in units of 10 6 g cm -3 
I a = a function of p, Z, A. 



(29) 
(30) 
(31) 



The temperature or density of the cross-over from phonon 
dominated to impurity dominated process depends on the 
impurity strength Q, given by, 



= ~y^ni(Z-Zi) 2 
n ^ 



(32) 



where n is the total ion density, n, is the density of impu- 
rity species i with charge Zi, and Z is the ionic charge in 
the pure lattice (Yakovlev & Urpin 1980). 

It should be mentioned here that in an accretion heated 
crust the effect of the impurities is negligible except when 
the rate of mass accretion is small implying a lower crustal 
temperature. This is true even when the crustal currents are 
concentrated at higher densities. To demonstrate this, we 
present some calculations for different values of Q in fig.(|l|). 
In this plot the evolution of the surface field with time has 
been shown for two different values of the impurity con- 
centration. The calculations pertain to two different rates 
of accretion and for currents concentrated at two different 
densities. It is evident from the curves that the impurity 
content is important only for lower rates of mass accretion. 
It should be noted that in fig.([j]) the curves 2Aa, 2 Ah and 
2Ba, 2Bb are actually pair-wise indistinguishable signify- 
ing that the effect of the impurity content is completely in- 
significant for somewhat larger rates of accretion. Since, in 
the present work we concentrate only on demonstrating the 
effect of space-time curvature on the evolution of the mag- 
netic field, Q = is assumed for the rest of our calculations. 
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Figure 2. The evolution of the surface magnetic field in an ac- 
creting neutron star of mass 1.4 Mq corresponding to an interme- 
diate equation of state given by WFF. Curves 1 to 3 correspond 
to densities of 10 13 , 10 12 , 10 11 gem -3 respectively, at which the 
initial current configurations are centred, assuming a flat space- 
time. Whereas, curves la to 3a correspond to curved space-time. 
All curves correspond to M = 10 -10 M0/yr and Q = 0. 



Figure 3. The evolution of the surface magnetic field in an ac- 
creting neutron star of mass 1.4 MQCorresponding to an interme- 
diate equation of state given by WFF. Curves 1 to 3 correspond 
to M = 10~ 12 , 10" 11 , 10~ 10 Mo/yr respectively, assuming a flat 
space-time. Curves la to 3a correspond to curved space-time. All 
curves correspond to Q = and a density of 10 11 gem -3 at which 
the current distribution is initially concentrated. 



5 RESULTS AND DISCUSSIONS 
5.1 Final Surface Fields 

We summarise the results of our investigation in the follow- 
ing plots. In fig.(Q) we plot the evolution of the surface mag- 
netic field for various densities of current concentration and 
for a fixed rate of mass accretion. It is clearly seen that the 
difference in final field values between that obtained using a 
flat and a curved space-time is ~ 1/4 orders of magnitude 
whereas that obtained by using two different densities for 
the initial current concentration could be as large as an or- 
der of magnitude. Similarly in fig.(j3|) the difference in the 
final field values obtained using two different rates of accre- 
tion is much larger than any modification introduced by a 
general relativistic framework. Since the values of the cur- 
rent concentration density or the rate of accretion used by 
us are well within the limits of observational uncertainty the 
inclusion/exclusion of a general relativistic framework does 
not really affect the calculation of the evolution of crustal 
magnetic field. 

In fig.s (Q) and (^[) we compare the effect of curved space- 
time for neutron stars of different masses and for different 
equations of state. The actual amount of field decay is more 
for a compact (a stiffer equation of state or a larger mass) 
star owing to a thinner crust and hence a smaller diffusive 
length-scale. But the modification due to general relativistic 
effects is larger for compact stars. That effect is discernible 
in fig.(^). But the effect of a mass difference is rather in- 
significant as is seen in fig.(^J). It should be noted here that 
the crustal models used in our calculation (described in sec- 
tion [0]) are somewhat older. The newer models predict 
lower crustal mass (Lorenz, Ravenhall & Pethick 1993) im- 
plying a higher surface red-shift for a given total mass of the 
star. It is evident from fig.(|^) that such an effect is again 
small compared to other factors mentioned above. 




log (t/yr) 

Figure 4. The evolution of the surface magnetic fields in neutron 
stars accreting at a rate M = 10 -10 Mo/yr, corresponding to an 
intermediate equation of state given by WFF. Curves 1 and 2 
correspond to the total masses of the neutron star being equal to 
1.4 and 1.8 Morespectively, assuming a flat space-time. Curves 
la and 2a correspond to curved space-time. All curves correspond 
to Q = and a density of 10 11 gem -3 at which the current 
distribution is initially concentrated. 



5.2 Effect of the Curvature of Space-Time 

To understand the effect of the curvature of space-time on 
the evolution of a purely crustal field we first define the char- 
acteristic diffusive and advective time-scales for eq.s jl4| ) 
and (H) as follows: 



diff 
Tflat 



diff' 
^"curved 



4naL 2 



- 1 ' 2 47RTL 2 



adv 
''"flat. 



adv 
''"curved 



L 

v(r) ' 

L 



,V2 



(33) 
(34) 
(35) 
(36) 
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log (t/yr) 

Figure 5. The evolution of the surface magnetic fields in neutron 
stars, of mass 1.4 Mq, accreting at a rate M = 10 -10 Mo/yr. 
Curves 1 and 2 correspond to W74 and P71 equations of state 
respectively, assuming a fiat space-time. Curves la and 2a cor- 
respond to curved space-time. All curves correspond to Q = 
and a density of 10 11 gem -3 at which the current distribution is 
initially concentrated. 



M/Mq 


P71 


WFF 


W74 


1.4 


1.499 


1.263 


1.23 


1.6 




1.322 


1.278 


1.8 




1.396 


1.335 


2.0 






1.408 



to the masses and the equations of state used in table 



cspc 



where L is the length-scale associated with the underlying 
current distribution supporting the field. The extra factor in 
T dlff for the curved space-time is, of course, the well known 
red-shift factor introducing an overall slow down compared 
to the case of a flat space-time. In table (|^) we compare the 
values of the red-shift at the surfaces of the stars of differ- 
ent masses and for different equations of state. This factor 
is larger for a stiffer equation of state and for a larger mass 
- both of which makes the star more compact. But the nu- 
merical value of the red-shift factor is not very different from 
unity as can be seen from table (Q) . It should be noted that 
in defining r adv we have simply collected terms proportional 
to the first space derivatives of g(r,t). Hence, the extra term 
in the numerator of r adv for the curved space-time is actu- 
ally a modification to the diffusive part due to the curvature 
of space-time and not really a modification to the convec- 
tive part. But this definition allows us to make upper bound 
estimates for the modifications to the final field values, due 
to the curvature of space-time, as described below. 



It is evident from the above definitions that: 

diff - diff i adv . adv 

r flat ^curved ana Tfl at > T cm . vcd . 



(37) 



This indicates that while the red-shift factor introduces an 
overall slow-down in the field evolution, the extra term, pro- 
portional to d r g(r,t) arising out of the curvature of space- 
time, accelerates the decay. Therefore, the actual modifica- 
tion in the final surface field is given by the competition of 



these two effects. In order to make an estimate of this dif- 
ference between the final surface field values obtained using 
a flat and a curved space-time, let us assume that g(r,i) 
evolves in time with an approximate exponential behaviour 
with the above time-scales. In that case, the difference is 
obtained as: 



A = 



, -Bfflflat 
'B(0)flat 



- log 



B(t) c 



S(0)c 



, -BMflat 
'S(0)flat 



= 1 



Tflat 
^"curved 



(38) 



for either diffusion or convection. Before making estimates 
using the above definition we should note that the time evo- 
lution of g(r,t) does not really follow an exponential be- 
haviour, particularly at late times. Therefore, the estimates 
that we make using these time-scales are nothing but very 
crude upper limits. 



For a star of mass 1.4 M0(with an eos given by WFF, 
with the currents concentrated initially at a density of 
10 11 gcm~ 3 and accreting at a rate of 10 -10 Mq yr _1 we 
find that Adiff ~ 20% and A a dv ~ 5%. This implies that, at 
most, we should see a difference of ~ 15%. Consider curve 
3 and 3a in fig.(Q) (which corresponds to these parameters). 
We see a difference of about half an order of magnitude (be- 
fore the curves level off and the effect of inner boundary 
condition of 'flux-freezing' starts controlling the evolution) 
which is close to ~ 10% of the corresponding value of the 
surface field and therefore within the upper limits estimated 
above. The actual difference between the final surface fields 
are, of course, much smaller and is closer to ~ 5%. 



6 CONCLUSIONS 

From the results presented above we conclude that the effect 
of the curvature of space-time on the evolution of a crustal 
magnetic field of accreting neutron stars is quite small. To 
summarise: 

(i) the difference in the final field strength obtained by 
using a general relativistic calculation is much smaller com- 
pared to the uncertainties already existing due to the plausi- 
ble range of various crustal and binary evolution parameters; 

(ii) the effect of curvature of space-time is larger for a 
stiffer equation of state or for a massive star but such dif- 
ferences are again quite insignificant. 
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